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^ '. Abstract 
^ . 

' We construct static solutions to a SU{2) Yang-Mills (YM) dilaton model in 4 + 1 

. dimensions subject to bi-azimuthal symmetry. The YM sector of the model consists 

I I of the usual YM term and the next higher order term of the YM hierarchy, which 

Qj ' is required by the scaling condition for the existence of finite energy solutions. The 

■ basic features of two different types of configurations are studied, corresponding to 
^ , (multi)solitons with topological charge n , and soliton-antisoliton pairs with zero 

topological charge. 
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1 Introduction 

Multi-instantons and composite instanton-antiinstanton bound states subject to bi-azimuthal 
symmetry were reported in a recent paper |1] . These were constructed numerically, for the 
usual {p = 1) SU{2) Yang-Mills (YM) system in 4 Euclidean dimensions, the spherically 
symmetric special case being the usual BPST [2] instanton. 

In a work [3j unrelated to [Ij, regular and black hole static and spherically symmetric 
solutions to a Einstein-YM (EYM) system in 4-1-1 dimensional spacetime were constructed 
numerically. The YM system in that model had gauge group 50(4), with the connection 
taking its values in (one of the two) chiral spinor representations of 5*0(4), namely in 
SU{2). Given that the solutions in [3] were static, i.e. that the YM field is defined on a 
4 dimensional Euclidean space, the SU{2) YM field in [3j is the same one as that in [1]. 
This is the relation between the two works [T| and [3J, and our intention here is to exploit 
this relation. 
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The present work serves two distinct purposes. The first and main purpose is to pave 
the way for the construction of more general, non-spherically symmetric solutions to EYM 
systems in a five dimensional spacetime. To our knowledge, no such results in EYM 
theory have appeared in the literature to date. Although considerable progress has been 
made in constructing asymptotically flat higher dimensional EYM solution^ all known 
configurations were subject to spherical symmetry. Our choice of a YM-dilaton (YMd) 
model is made because it has been shown that, at least in = 3 + 1 dimensions, the classical 
solutions of this system mimic the corresponding EYM solutions so the dilaton-YM 
exercise serves as a warmup for the considerably more complex gravitational problem. 
Our choice of a YMd model is made as an expedient in attempting the analysis of the 
corresponding EYM model, the last being of physical interest low energy effective actions 
of string theory, descended from 11 dimensional supergravity. It is also a coincidence 
here, that these supergravity descended low energy effective actions include the dilaton 
in addition to gravities and non Abelian matter. But here, the dilaton appears only as a 
substitute for gravity. 

A particular feature of the model to be introduced in the next section is that it features a 
term that is higher order in the YM curvature. As will be explained in section 2, such terms 
are necessary to ensure that the solution yields a finite mass. Such terms were employed in 
previous works [H IH [5l El [7l E] with precisely the same purpose. The physical justification 
for introducing higher order YM terms, which goes hand in hand with the inclusion of 
higher order gravitational terms, is that these occur in the low energy effective action of 
string theory [9] . Thus in principle the choice of higher dimensional EYM models involves 
the selection of higher order terms in the gravitational and non Abelian curvatures, namely 
the Riemann and the YM curvatures, which are reparametrisation and gauge invariant. 
Because we are concerned with finding classical solutions, we impose a pragmatic but 
important further restriction, namely that we consider only those Lagrange densities that 
are constructed from antisymmetrised 2p curvature forms, and exclude all other powers 
of both Riemann and YM curvature 2-forms. (In the gravitational case this results in 
the familiar Gauss-Bonnet type Lagrangians, while in the case of non Abelian matter to 
the YM hierarchy pointed out in footnote 3 below.) As a result, only velocity-squared 
fields appear in the Lagrangian, which is what is needed both for physical reasons and for 
solving the classical field equations. In practice we add only the minimal number of such 
higher order terms that are necessitated by the requirements of finite mass. This criterion 
makes the inclusion of higher order gravitational terms unnecessary since we know from 
the (numerical) results of [1] that the qualitative properties of the classical solutions are 
insensitive to them. In addition to this argument based on numerical results there is an 
independent argument advanced at the end of section 2 of [6], based on the symmetries of 
the (higher order) gravitational terms, which in the absence of a dilaton dispenses with the 

^ EYM particle-like and black hole solutions approaching at infinity the Minkowski background have 
been constructed numerically for d = 6,7 and 8 [4]. The properties of globally regular solutions in arbitrary 
dimensions have been studied both numerically and analytically in |^5j . Higher dimensional asymptotically 
(anti-)de Sitter solutions have been found numerically in [6], [7], as well as such systems whose gravitational 
sector consists of higher order Gauss-Bonnet like gravitational terms [8] . 
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effectiveness of employing such terms. (Note here that the present YMd model is being 
used as a prototype for a EYM model, without an additional dilaton field.) This leaves 
one with higher order YM curvature terms only, whose status in the context of the string 
theory effective action is complex and as yet not fully resolved. While YM terms up to 
arise from (the non Abelian version of) the Born-Infeld action [10], it appears that this 
approach does not yield all the terms [11]. Terms of order F^ and higher can also be 
obtained by employing the constraints of (maximal) supersymmetry p2]. The results of 
the various approaches are not identical. In this background, we restrict our considerations 
to terms in the YM hierarchy (see footnote 3) only, in particular to the first two terms. 

Concerning our particular choice of 4 + 1 spacetime dimensions here, our reasons are: 
When imposing axial symmetry on a YM field in d = D + 1 dimensions the simplest 
way is, following [T3], to impose spherical symmetry in the D — 1 dimensional subspace 
of the d spacelike dimensions. In this case the Chern-Pontryagin topological charge is 
fixed by the boundary conditions imposed on the first polar angle, and no analogue of the 
vortex number appearing in the axially symmetric Ansatz for d = 3 [12] is featured [2]. 
Technically, the absence of a vortex number makes the numerical integration much harder. 
Imposing axial symmetry in turns in the x — y and z — u planes of D = 4 Euclidean space 
as in [1] on the other hand, features two (equal) vortex numbers, making the numerical 
work technically more accessible. It is our intention to use the particular bi-azimuthally 
symmetric Ansatz of pj in D = 4 that has led us to restrict ourselves to = 4 + 1 
dimensional spacetime. (Numerical work on implementing axial symmetry like in [H] is 
at present in active progress.) Of course, the exploitation of this type of symmetry is not 
restricted to 4 + 1 spacetime, but can be extended to any odd 2q + 1 spacetime where q 
distinct azimuthal symmetries are imposed, but this in practice results in residual PDE's 
of order three and higher for g > 3. 

The second and subsidiary aim of this work is to break the scale invariance of the usual 
YM system in D = 4 studied in [Ij, and the introduction of the dilaton field does just 
that. The question of instanton-antiinstanton bound states in a scale breaking model is an 
interesing enough matter in itself, presenting a second important motivation for this work. 

In Section 2 we present the model, impose the symmetry and state the boundary con- 
ditions, in successive subsections. The numerical results are presented in Section 3, pre- 
senting both solutions with spherical and bi-azimuthal symmetry. We give our conclusions 
and remarks in the final section. 



2 The model 

The model in 5 spacetime dimensions with coordinates xm = {xo,Xfj^) that we study here 
is described by the Lagrangian 

2 I / ''"l 2a4, -772 I "^2 6a0 rp ^2 
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where is the dilaton field, J^mn = f^M-^Af — f^Jv-^A/ + [-4A/7-^Af] is the 2-form YM cur- 
vature and J^MNRS = {^M[N, ^RS]} is the 4-form YM curvature consisting of the totally 
antisymmetrised product of two YM 2-form YM curvatures. (The bracket [i^pcr] implies 
cyclic symmetry.) ti and T2 are dimensionful coupling strengths which will eventually be 
scaled out against the constant a in the exponent, which has the inverse dimension of the 
dilaton field 0. Similar to the d = 3 + 1 case, the form we choose for the coupling of the 
dilaton field to the nonabelian matter was found by requiring that a shift ^ + 00 of 
the dilaton field to be compensated by a suitable rescaling of the coordinates. 

Let us give a brief justification for the choice of the model ([T]). At the most basic 
level it is a YM-dilaton (YMd) model designed to simulate qualitatively a EYM model in 
d = 5. Repacing the dilaton in ([1]) by a gravitational term is a physically relevant model, 
representing part of a low energy effective action in d = 5. Adding gravitational terms to 
([1]) as it stands is a EYMd model, which is just as physically relevant. 

The YM system, which scales as L~^, in c? = 4 + 1 supports static solitons, namely 
the BPST instantons in d = 4 + dimensions. When the usual Einstein-Hilbert gravity, 
which scales as L~^, is added to the YM term, the soliton collapses because of the (Derrick) 
scaling mismatch. To compensate for this scaling mismatch, a term scaling as L~'^, with 
z/ > 5 must be added. If one is to restrict to positive definite terms 0, u will be even, and 
the most economical choice is u = 6. A typical such term would be Tr {F A DX)^, where 
X is a scalar field, e.g. a Higgs or sigma-model field. This necessitates the introduction 
of a completely new (scalar) field which unlike the dilaton is not directly recognised as a 
constituent of a low energy effective action. For this reason we eschew this choice, and 
restrict our attention instead to systems featuring only YM (and eventually YMd) fields. 
The most economical choice then is to compensate with the YM term Tr (F A F)^, scaling 
with u = 4:. We note, finally, that adding a (positive or negative) cosmological constant 
does not remedy the scaling mismatch since these terms do not scale at all. Indeed, in all 
higher dimensional EYM cases studied, with A = 0[il|3|,A<0[6j and A > [7J, the 
mass turns out to be infinite when higher order YM terms are not employed. 

The YM sector of the action density in 4 -(- 1 dimensions employed here, is that one 
used in [3] , namely the superposed p = 1 and p = 2 members of the YM hierarch}{§. The 
dilaton breaks the scale invariance of the usual p = 1 YM system, and a simple Derrick- 
type scaling argument shows that no finite mass/energy solution can exist if T2 = 0, i.e. 
the p = 2 term in ([1]) is necessary. 

The YM and dilaton field equations read 



^The choice of a Chern-Simons term e\f_i^paTrA\ [Ff^„Fpcr ~ F^^ApA^ + ^ApAi,ApAcr) scaling as L^^ is 
a possibility, albeit a considerably harder problem technically, and is at present under active consideration. 

■^The YM hierarchy labeled by the integer p was introduced in [TB] in the context of self-dual solutions 
in Ap Euclidean dimensions, but superpositions of various p members were employed ubiquitously since. 



(2) 




(3) 



4 



In ([3]) we have used the notation 



-^1 — ^~;^Tr:F|^^ ' -^2 — ^^^Tr:F^^^^5 . (4) 
2.1 Imposition of symmetry and residual action 

In the YM connection Am = {-Ao, A^), we choose the temporal component = to vanish 
and the spacelike components A^ is subjected to two successive axial symmetries, described 
in [1]. We denote the Euclidean four dimensional coordinates as = (x, y; z, u) = {xa, Xt), 
with a = 1,2 and i = 3, 4, and use the following parametrisation 

Xa = r sin 9 Xa = pxa, Xi = r COS 9 Xi ^ axi , (5) 

where = = |xq,P + with the unit vectors appearing in ([5]) parametrised as 
Xa = (cosyji, sin(y9i), Xi = (cos (/?2, siu (/?2) , with < 6* < I spanning the quarter plane, 
and the two azimuthal angles < y^i < 27r and < < 27r. The d = A + 1 Minkowski 
spacetime metric reads for this ansatz 

ds^ = -de + dr^ + r^{d9^ + sin^ 9d^l + cos^ 9d^l). (6) 

The first stage of symmetry imposition is of cylindrical symmetry in the x^ = {xi,X2) 
plane, and the Ansatz is stated as is 

A = (^) + (^) (£n(i))^ + Af - Af xjl^ , (7) 

A = AT'>nfT.p^-AfT.u , (8) 

in which the index m = 3, 4 is summed over, and the unit vector ^ = (cos^iy^i , sin niipi) 
is labeled by the vortex integer rii, Eap being the Levi-Civita symbol. The spin matrices 
E^i, = Eqj, Ejj) in (171) and (IHD are one or other of the two chiral representations of 

5*0(4), i.e. they are SU{2) matrices. 

If in ©-(IE]) we regard the functions (0™, 0^) = 0" as a 5*0(3) isovector field, and 
{A'^^^Af) = Af and {A'^^^Af) = Af as 50(3) YM connection fields, with antisymmet- 
ric 50(3) algebra indices [ab] = ([34], [45], [53]), then it turns out that J-'^^ = {J-'a/s, ^ai, ^ij) 
is expressed exclusively in terms of the curvature {F°--',F°-p) of the 50(3) connection 
{A'^p^Af'), and the corresponding covariant derivative {Di(j)°', Dp(j)°') of 0", all defined on 
the hyperbohc space with coordinates {xi,Xp). 

The second stage of symmetry imposition is expressed most succintly by rewriting the 
residual fields (A^'', A^'^) and 0" on this hyperbolic space in matrix representation 

A = ~Af^,, , Ap = -^-Af^,,, (9) 

$ = r^a4. (10) 

Azimutal symmetry in the (xs — X4) plane is imposed on the connection fields by 
decomposing At formally according to ^^^| and Ap according to ([8]). Noticing now that the 
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index a in (p!)- (fTOl) runs only over three values, and reassigning the values of the index 
i = 1,2, the analogues of ([7]) and ([8]) contract to give 

A = (^^^^^^ j:^,,x, + (^^^ (11) 



A, = Al^nfE,^, (12) 

exhibiting the Abelian connection A^^ analogous to the Abelian connection A^^ appearing 
in ([7j) and the isodoublet function (x^;X^)- 

The corresponding axially symmetric decomposition of $ in fllOj) is 

$ = ^1^^ S,4 + e'S34. (13) 

In (fTT]) - (IT^ and (IT^ we have used the unit vector nf'^ = (cosn2V52 , sin 712(^2), with vorticity 
integer n2. The final stage of symmetry imposition is to treat the two azimuthal symmetries 
imposed in the x — y and the z — u planes on the same footing, leading to the equality of 
the two vortex numbers, rii = n2 = n. 

Denoting the residual functions {A^p^,A^^) = {apa^,), (x^^x'') = X^, (^\^^) = 
and regarding ( Abelian connection on the quater plane defined by (p, a) , the 

residual action densities can be expressed exclusively in terms of the SO (2) curvature 

fpa = dpa„ -d„ap (14) 

and the covariant derivatives 

^^pX^ = '9pX^ + «p(^x)^ , P,X^ = 9.X^ + a,(£X)^ (15) 

The residual two dimensional YM action densities descending from the p = 1 and the p = 2 
terms Li and L2 defined by (jlj) are, respectively. 



cr ' p ^ 'pa 



12 pa 



[sABX^e f,. + T)^,X^T)^\e) . (16) 



These residual action densities are scalars with respect to the local 5*0(2) indices A,B, 
hence they are manifestly gauge invariant. They describes a U{1) Higgs like model with 
two effective Higgs fields ^"^^ ^"^1 coupled minimally to the U{1) gauge connection 
(ap,ao-). To remove this U{1) gauge freedom we impose the usual gauge condition 

dpUp + d^a^ = . (17) 
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Since our numerical constructions will be carried out using the coordinates (r, 9) we display 
(fT6|) also as 
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rsinecosej%+'-^(\V^Xr + \\VeX^\' 
cos o \ 

rcosO ( ,^ ^4,9 1 ,^ ^A^'\ 



1 



sin0 \ / rsin^cos^ 



12 sin Q cos ^ 
The total mass-energy M of the system is 



(18) 
(19) 



M 



oo /"''"/2 

m — I dr dO 
io 



^r=^sin^cos^(02^ + ^0^^) + (e^^^^Li + e^^-^Ls) 



,(20) 



and equals the total action of solutions, viewed as solitons in a c? = 4 Euclidean space. 
2.2 Boundary conditions 

To obtain regular solutions with finite energy density we impose at the origin (r = 0) the 
boundary conditions 



Or = , ae = , x' 





-n2 





-rii 



(21) 



which are requested by the analyticity of the YM ansatz, and c}r0|r=o = for the dilaton 
field. In order to find finite mass solutions, we impose at infinity 



Oj. = 0, 



-2m, x' 



1)^+1^2 



sin 2m9 
cos 2m9 



sin 2m9 
cos 2m9 



0, 

(22) 



m being a positive integer. Similar considerations lead to the following boundary conditions 
on the p and a axes: 

ar = —dri\ ae = —dee, = 0, = 0, = 0, = -m, 9^0 = 0, (23) 

rii rii 

for 6^ = 0, and 

a, = -9,x\ ae = -9ex\ x' = 0, = 0, = -^2, 9^^' = 0, 9,0 = 0, (24) 
n2 n2 



for 9 = 71/2, respectively. 
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2.3 Topological charge 



In our normalisation, the topological charge is defined as 



Q — 22^2^'^'^'^°' / '^^{•^ ui^-^ o(t\ d X, (25) 



which after integration of the azimuthal angles {ipi,ip2) reduces to 

<1 = \^^^ [ d^ABX^^'' U. + 'D,x''V.^A d'x (26) 



2 / V2 



1 1 e,,d,{x^V,e-i^'D,X^)d^x. 



(27) 



The integration in fl26l) is carried out over the 2 dimensional ). As 

expected this is a total divergence expressed by fl27j) . 

Using Stokes' theorem, the two dimensional integral of fl27j) reduces to the one dimen- 
sional line integral 

q = \ I X^V.eds,, (28) 



This integral has been evaluated in [T] by reading off the appropriate values of x"^ and 
from ([23])-([2lD. The result is 

q=\[l- {-ir]n^n,. (29) 



3 Numerical results 

Apart from the coupling constants Ti and T2 the model contains also the dilaton constant 
a. Dimensionless quantities are obtained by rescaling 

r-^r{T2/T,Y'\ (30) 

This reveals the existence of one fundamental parameter which gives the strength of the 
dilaton-nonabelian interaction 

2 2 3/2/ 1/2 /01 \ 

a = a T^' jr^ , (31) 

which is a feature present also in the EYM case j3]. We use this rescaling to set ri = 1, 
r2 = 1/3 in the numerical computation, without any loss of generality. 

One can see that the limit a — > can be approached in two ways and two different 
branches of solutions may exist. The first limit corresponds to a pure p = 1 YM theory with 
vanishing dilaton and p = 2 YM terms, the solutions here replicating the (multi-)instantons 
and composite instanton-antiinstanton bound states discussed in pQ. The other possibility 
corresponds to a finite value of the dilaton coupling a as ri — 0. Thus, the second limiting 
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configuration is a solution of the truncated p = 2 YM system interacting witli tlie dilaton, 
witli no p = 1 YM term. 

We liave studied YMd solutions with m = 1, 2. From our knowledge of the tolopogical 
charges (|29i) . the m = 1 solutions will describe (multi)solitons and the m = 2 solutions, 
soliton-antisoliton configurations. Also, to simphfy the general picture we set rii = n2 = n 
in the boundary conditions (l2T|) - flM|) . 

The spherically symmetric solutions are found by using a standard differential equations 
solver. The numerical calculations in the bi-azimuthally symmetric case were performed 
with the software package CADSOL, based on the Newton- Raphson method In this 
case, the field equations are first discretized on a nonequidistant grid and the resulting 
system is solved iteratively until convergence is achieved. In this scheme, a new radial 
variable x = r/(l + r) is introduced which maps the semi-infinite region [0, oo) to the 
closed region [0, 1]. 

As will be described below, solutions exist for certain ranges of the parameter a. It 
turns out that m = 1 solutions with all n and m = 2 solutions with n = 1 exist for a 
range of a starting from a a — > limit, but do not persist all the way up to the second 
a limit. (However, the way the solutions approach the limit a —>■ depends on m.) 
By contrast we find that m = 2 solutions with all n > 1, exist for all a between the two 
limits. 



3.1 m = 1 configurations 

n = 1 spherically symmetric solutions 

In the spherically symmetric limit, which case we shall analyse numerically first, the angular 
dependence of these functions is fixed and the only remaining independent function depends 
on the variable r. The independent function in this case is ae = w{r) — l, with the remaining 
fuctions {aryX^y^"^) given by 

= 0, = = -(ti;(r)-l)sin2^, = -(w(r)-l) cos^ ^-1, = -(u;(r)-l) sin^ ^-1. 

^ (32) 
The functions (/)(r) and w{r) satisfy the equations 

\ r r-^ ) 



/2 



e'^'rw'in + 3e^'^r2 ^^ , ^ ) = ^ n + 3r2e^^^ . (33) 

The asymptotic solutions to these functions can be systematically constructed in both 
regions, near the origin and for r ^ 1. The small r expansion is 

w{r) = l- br^ + 0{r^), = 0o + 4a^( y + 9r2&^)6V + 0(0, (34) 
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with b, 00 two real parameters, while as r — oo we find 



w{r) = ±1 



401 



40? 



+ 0(l/r«), 



01 32aV^e2'^o 



+ 0(l/O. 



(35) 



27r6 



27r6 



We numerically integrate the Eqs. ( l33i) with the above set of boundary conditions for 
Ti = 1,T2 = 1/3 and varying a. The picture we found is very similar to that found for the 
EYM system [3], the dilaton coupling constant playing the role of the Newton constant. 
First, for a given a, solutions with the right asymptotics exist for a single value of the 
"shooting" parameter b which enters the expansion (IMI) . For a small enough, a branch of 
solutions smoothly emerges from the BPST configuration [2] . When a increases, the mass 
M and the absolute value of the dilaton function at the origin increase, as indicated in 
Figure 1. These solutions exist up to a maximal value amax — 0.36928 of the parameter a. 

As in the corresponding gravitating case [3], we found another branch of solutions in 
the interval a G [acr{i),Cimax] with a^^^^-^ ~ 0.2653. On this second branch of solutions, 
both 0(0) and M continue to increase but stay finite. However, a third branch of solutions 
exists for a G [0.2653,0.2652], on which the two quantities increase further. A fourth 
branch of solutions has also been found, with a corresponding acr(3) — 0.2642. The mass 
M, the value of the dilaton field at the origin 0(0) and the initial (shooting) parameter b 
increase along these branches. Further branches of solutions, exhibiting more oscillations 
around a ~ 0.264 are very likely to exist but their study is a difficult numerical problem. 
This critical behaviour is described as a conical fixed point in the analytic analysis in |5] . 
Therefore we conclude that, as in the spherically symmetric gravitating case [3], the limit 
Ti = is not approached for solutions with m = 1, n = 1. 

As a general feature, all solutions discussed here present only one node in the gauge 
function w{r). As in the higher dimensional EYM models discussed in [HIS], no multinode 
solutions were found. 



Solutions with bi-azimuthal symmetry with nontrivial dependence on both r and 6 are 
found for {ni,n2) ^ 1 subject to the boundary conditions fl?Il) - fl2^ . We have studied 
solutions for m = 1 with 2 < n < 5. The general features of the m = 1 solutions are the 
same for all n > 1. Also, as seen in (129|) . the m = 1 configurations carry a topological 
charge q = n^. The corresponding solutions of the J-'^j^ model are self-dual and have been 
considered already in [12], [IS] (for a different parametrization of the gauge field, however). 

These solutions are constructed by starting with the known spherically symmetric con- 
figuration and increasing the winding number n in small steps. The iterations converge, 
and repeating the procedure one obtains in this way solutions for arbitrary n. The physical 
values of n are integers. The typical numerical error for the functions is estimated to be 
of the order of 10~^ or lower. 

Any spherically symmetric configuration appears to result in generalisations with higher 
winding numbers n. Moreover, the branch structure noticed for the m = 1, n = 1 case 
seems to be retained by the higher winding number m = 1 solutions. Again, the first branch 



n> I 
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of solutions exists up to a maximal value of a, where another branch emerges, extending 
backwards in a. We managed to construct higher winding number n counterparts of the 
first two branches of spherically symmetric solutions. The mass M and the absolute value 
of the dilaton function at the origin increase along these branches, as shown in Figure 2. 
Note that the value of the dilaton function at the origin exhibited in the figures is actually 
0(r = 0, ^ = 0), restricting to = 0. This restriction is reasonable since for all solutions 
with bi-azimuthal symmetry discussed in this paper, the dilaton function at r = presents 
almost no dependence on the angle 6. 

We expect that the oscillatory pattern of 0(0) arising from the conical fixed point 
observed above for the spherically symmetric n = 1 solutions, will also be discovered for 
the n > 1 solutions here, but their construction is a difficult numerical problem beyond 
the scope of the present work. 

In Figure 3 we present the gauge functions, the dilaton, and the topological charge 
density 

read off fl26p . as functions of the radial coordinate r for five different angles for a typical 
first branch m = 1, n = 3 solution with a = 0.21. The functions ag and (f) have a small 
6 dependence (although the angular dependence increases with n), while xi ^^id .^i have 
rather similar shapes. The action density C possesses one maximum on the 6 = tt/A axis. 
All multicharge solutions found have concentrated energy and charge density profiles where 
individual (unit) charge consituents do not appear as distinct components. The moduli of 
the effective Higgs fields |x| = ix^X'^Y^'^ and |^| = {^'^^^Y^'^ possess one node each on the 
p and the a axes, respectively, which coincide with the maximum of the action density. 
The position of this node moves inward along the first and second branches. 

3.2 m = 2 configurations 

The m = 2 configurations can be thought of as composite systems consisting of two com- 
ponents which are pseudoparticles of topological charges ±n. Thus, these configurations 
reside in the topologically trivial sector and carry no Chern-Pontryagin topological charge. 
This type of solutions have no spherically symmetric limit. Also, their behaviour as a 
function of a is different from those with m = 1 presented above, in that solutions for all 
values of a exist between the two distinct limits of a — implied by (13T|) . for all n except 
for n = 1. 

n = 1 

It is perhaps interesting to note from the outset that m = 2 , n = 1 solutions to be 
described now, have apparently no counterpart in the 4 + dimensional p = 1 YM model 
studied in |1]. (It turns out that for the (m = 2,?t, = 1) solution in that case there is no 
analytic proof of existence either.) The obvious difference of the 4 + 1 dimensional model 
([1]) here and the 4 + dimensional p = 1 YM model is that the solutions of the former are 
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parametrised by the effective coupling constant a, while the latter has no such parameter. 
As will be described below, m = 2 n = 1 solutions exist for a certain range of a, and this 
range excludes the limiting case where the contribution to the action of the dilaton term 
and the p = 2 YM term in ([T]) disappear, i.e. a J-'lj^ model. 

We find that in the limit a — >■ resulting from a —>■ 0, cf. (1311) . no solutions of this 
type exist. However in the limit a ^ corresponding to a finite value of the dilaton 
coupling a as Ti ^ 0, such solutions exist. This limiting configuration is then a solution 
of the truncated system consisting of the dilaton term and p = 2 YM term J^mj^jis^ which 
dominate. Its characteristic feature is that for this configuration both nodes of the effective 
Higgs fields |x| and |^| merge on the 6 = tt/A axis. A family of solutions of the model ([T]) 
emerges from this configuration. As a increases, the nodes move towards the symmetry 
axes, p and a, respectively, forming two identical vortex rings whose radii slowly decrease 
while the separation of both rings from the origin increase. At the critical value ctcr — 0.265 
the node structure of the configuration changes, both vortex rings shrink to zero size and 
two isolated nodes appear on each symmetry axis. This structure is known for the usual 
YM system in d = 4 + [Ij, indeed, increasing of a along this branch can be associated 
with increasing of the coupling ti w.r.t. T2 as the dilaton coupling a remains fixed; then the 
term becomes leading. The maximum of the action density however is still located 

on ^ = 7r/4 axis. 

Another similarity with the instanton-antiinstanton solution of the d = 4 + p = 1 YM 
theory is that the gauge functions a^, ctg as well the dilaton function of the n = 1, m = 2 
solutions also are almost ^-independent. Along this branch the mass of the solutions grows 
with increasing a since with increasing coupling ri the contribution of the term J-'lj^^ also 
increases. 

As the effective coupling increases further beyond acr the relative distance between the 
nodes increases, one lump moving towards the origin while the other one moves in the 
opposite direction. Along this branch both the value of the dilaton field at the origin |0(O)| 
and mass of configuration M increase as a increases. This branch extends up to a maximal 
value amax — 0.311 beyond which the dilaton coupling becomes too strong for the static 
configuration to persist. The second branch, whose energy is higher, extends backwards 
up to amax — 0.279. Along this branch both |0(O)| and the mass of the configuration 
continue to increase as a decreases. Also the separation between the nodes decreases and 
both nodes invert direction of the motion, moving toward each other along this branch. 
In Figure 5 we present the values of the dilaton function at the origin 0(0) and the total 
mass (rescaled by ct^) of these configurations as functions of a. 

n = 2 

This configuration also resides in the topologically trivial sector and can be considered as 
consisting of two solitons of charges n = ±2. Then the interaction between the nonabelian 
matter fields becomes stronger than in the case of unit charge constituents and the expected 
pattern of possible branches of solutions is different from the n = 1 case above. 

Indeed, the n = 2 , m = 2 solutions show a different dependence on the coupling con- 
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stant a, with two branches of solutions. The lower branch emerges from the corresponding 
solution in pure p = 1 YM theory with vanishing dilaton and p = 2 YM terms, replicating 
the corresponding solution in [l]. The variation of the effective coupling along this branch 
is associated with the decrease of ri, at fixed T2 and fixed dilaton coupling a. The second 
branch emerges from a solution of the p = 2 YM-dilaton system, the unrescaled mass M 
diverging in this limit, with the rescaled mass Ma^ vanishing as seen from Figure 5a. At 
the maximal value amax — 0.2372 this branch bifurcates with the lower YM branch. For 
larger values of a, the dilaton coupling becomes too strong for the static configurations to 
persist. Thus for < a < Omax we notice the existence of (at least) two distinct solutions 
for the same value of coupling constant. 

For the same value of a, the mass of the second branch solution is larger that that of 
the corresponding lower branch configuration (s). One should also notice the existence of 
a curious backbending of the lower branch for 0.193 < a < 0.218. Four distinct solutions 
exist in this case for the same value of a (three of them located on the lower branch), 
distinguished by the value of the mass and the dilaton field at the origin. This pattern is 
illustrated in Figure 5. 

Again, observation of the positions and structure of the nodes of the effective scalar 
fields allows us to better understand the behaviour of the solutions. For lower branch 
solutions with small values of a there are two (double) nodes of the fields |x| and |^| on the 
p and a symmetry axes respectively. The locations of nodes correspond to the locations 
of the two individual constituents and the action density distribution posesses two distinct 
maxima on the 6 = 7t/4 axis. The distance between these nodes changes only slightly along 
the lower mass branch. The backbending in a observed in this case is reflected also for in 
the relative positions of the nodes. At the maximal value of a, the inner node is located 
at Po^'* = (Jo — 2.97 and the outer node is located at p'^^ = a^^ ~ 4.18. 

Along the upper branch, as a slightly decreases below amax, the inner node inverts 
direction of its movement toward the outer node which still moves inwards. Thus, both 
nodes on the symmetry axis rapidly approach each other and merge forming a two vortex 
ring solution as a ~ 0.2355. The action density then has a single maximum on ^ = 7r/4 
axis. As a decreases further both nodes move away from the symmetry axis and their 
positions do not coincide with the location of the maximum of the action density. Further 
decreasing a results in the increase of the radii of the two rings around the symmetry axis, 
and in the limit a — the rings touch each other on the 6 = tt/4 axis. 

In Figure 4 we give three dimensional plots of the modulus of the effective Higgs field 
^ for the n = m = 2 upper branch vortex solution at a = 0.20 and the n = m = 2 lower 
branch double node solution at the same value of a. The action density as given by ([T]) is 
also plotted at a = 0.20 both for the upper and for the lower branches. 

The numerical calculations indicate the possibility that the solutions of the fundamental 
YM branch, namely the branch on which the p = 1 YM term dominates, are not unique. 
It is possible that higher linking number configurations with higher masses might exist. 
This possibility will be explored elsewhere. 



13 



n = 3 



For the n = 3 configuration composed of two triple cliarged pseudoparticles we observe a 
somewliat simpler pattern. Tlie lower dilaton branch emerges from the limit of vanishing 
dilaton and p = 2 YM couplings and extends up to a maximal value acr = 0.165. Along 
this branch the configuration posesses two vortex rings. 

As a increases the mass of the solution increases and, at the same time, the radii of the 
rings slowly increase and both rings move inwards. This lower mass branch bifurcates at 
the critical value of the effective coupling with an upper branch which extends all the 
way back to a = (see Figure 5). Again, in this limit both vortex rings come into contact 
on the 6 = tt/4 axis. Thus, we observe no isolated nodes on the symmetry axis and both 
upper and lower energy branches correspond to vortex ring solutions. 

Note that the branch structure here closely resembles the pattern which was observed for 
the gravitating monopole-antimonopole chains and vortex solutions in d = 3 + 1 Einstein- 
Yang-Mills-Higgs system [201 [2T]. 

The profiles of a typical m = 2, n = 3 solution are presented in Figure 6 (the picture 
there applies as well for n = 1, 2 configurations). 

4 Summary and discussion 

Finite mass static solutions to a 4 + 1 dimensional SU{2) YMd model are constructed 
numerically. The YM sector of the model consists of the usual YM term, labeled p = 1, 
and the next higher order p = 2 term of the YM hierarchy [16]. The second YM term is 
necessary to counteract the scaling of the quadratic dilaton kinetic term, since we require 
finite energy solutions. In section 2 we have explained the rationale behind our choice of the 
model ([1]). Basically our criterion is that of compensating for the scaling of the quadratic 
kinetic dilaton term, by introducing an additional term that scales as , (z^ > 5), with 
the further criteria that this additional term be positive definite, and, that no new fields 
beyond the dilaton and the YM be employed. 

The solutions constructed are subject to the bi-azimuthal symmetry applied in [T]. 

Viewed as a 4 + dimensional model, this is a scale breaking version of a p = 1 YM 
theory whose bi-azimuthally symmetric solutions were presented in [1]. The latter model 
being scale invariant in d = 4 + 0, the present solutions have the new feature that they are 
referred to an absolute scale. This is an interesting feature of the present work. 

Our main motivation here, however, is to study a system which can give an insight 
into the qualitative features of static, finite mass solutions to a gravitating YM system in 
higher dimensions, which are not subject to spherical symmetry. To date the only higher 
dimensional Einstein- Yang-Mills (EYM) solutions known [H [3l [5l [6l [71 [8] are subject to 
spherical symmetry in the spacelike dimensions. Rather than tackle the appreciably more 
complex numerical problem of constructing non spherically symmetric EYM solutions, we 
consider here instead the corresponding YM-dilaton problem, knowing that the classical 
solutions of the latter simulate [13] the qualitative properties of the corresponding EYM 
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ones in (i = 3 + 1 dimensional spacetime. Thus the present work is a warmup to the ulti- 
mate aim of constructing non-spherically symmetric EYM solutions in higher dimensions. 
Having said that, we note that we have already made an appreciable start in the construc- 
tion of the corresponding d = 4 + 1 EYM solutions, and our results to date confirm the 
qualitative similarity of those with the dilaton-YM solutions presented here. Our results 
on the gravitating system will be reported elsewhere. 

In the context of giving a qualitative description of our results, it must be noted that 
in this first effort, we have restricted the dimensionality of the spacetime to 4 -|- 1. Thus 
we can only compare our results here with the 4-1-1 dimensional subset of the spherically 
symmetric EYM solutions [H El El El [H] , which were given in spacetime dimensions d > 5. 
We would expect however that the comparative features between spherically and the non 
spherically symmetric solutions which we uncover here, will stay qualitatively valid also in 
dimensions d > 5. 

As it happens, spherically symmetric EYM solutions in 4-1-1 dimensions [3j exhibit quite 
different qualitative properties compared to those in 5 -|- 1, 6-1-1 and 7-1-1 dimensions [1]. 
Indeed, the comparative patterns remain true modulo 4p dimensions, as explained in [5]. 
These features in question concern the branch structure of the said solutions with respect 
to the effective coupling parameter a in the problem. It turns out that for an appropriate 
EYM model in d = Ap + 1 dimensions, a peculiar branch structure, absent in d ^ Ap + 1, 
occurs. This was explained in [5] to be due to the occurence of what was called there a 
conical fixed point singularity, which manifests itself by the oscillatory behaviour in a of 
the global quantities near the critical value of the effective coupling parametei0. 

Now in the present work we do not employ gravity but have instead the dilaton, which 
is represented by the singlet scalar field 0. It is therefore unavoidable that the role of the 
metric function at the origin in the gravitating case, should be replaced here by 0(r = 0), 
the dilaton function at the origin. This correspondence can be made uniquely in the 
special case of our (m = 1 , n = 1) bi-azimutal solutions, which are simply the spherically 
symmetric subset. As seen from Figure 1, the oscillatory pattern of 0(0) is clear. 

Concerning (m > 2 , n) solutions, these are not spherically symmetric and depend on 6 
in addition to r. Instead we have employed the values 0(r = 0, ^ = 0) to track the branch 
behaviour of our solutions. This is quite a reasonable criterion, since we find that the 
dependence of 0(0, 6) is very small on all branches constructed. The branch structures for 
(m = 2 , n = 1, 2, 3) are displayed on Figure 5. We see that (m = 2 , n = 2, 3) solutions 
have largely similar patterns, except for the additional backbending showing up in the 
n = 2 case. The branch structure for (m = 2 , n = 1) on the other hand is drastically 
different. Indeed it seems quite reminiscent of the corresponding spherically symmetric 
(m = 1 , n = 1) solutions, pehaps exhibiting a conical fixed point behaviour too, but the 
numerical accuracy is entirely inadequate to decide this, either way. 

Our results concerning m > 2 solutions are here restricted to the concrete construction 

This feature persists also when a scalar matter field is added, for example when a gauged Grassmannian 
sigma model field is included in the Lagrangian [22]. We expect it to persists also when a Higgs field is 
added instead, in 4p + 1 dimensions. 
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of m = 2 solutions. However, our preliminary numerical results indicate that most quali- 
tative features remain true for m = 3, 4. Also, for m = 2, we have restricted to n = 1, 2, 3 
only, which is quite adequate. An interesting remark on m = 2 solutions is the existence 
of the (m = 2 , n = 1) confirmed by our results, which we had not found for the (1 = 4 + 
dimensional p = 1 YM model in [1] , and for which the analytic proof of existence is also 
absent. This is easy to understand since our (m = 2 , n = 1) solutions exist only for 
that range of the parameter a for which all terms in the Lagrangian ([1]) contribute to the 
action, and for the limiting range of a for which only the usual YM term {p = 1 term) 
would dominate, there are no solutions. So there is no special case of our solutions which 
could describe (m > 2 , n = 1) bi-azimuthal instantons that are absent in [1]. 

Apart from the above features, we have found a very rich pattern of the zeros of the 
moduli of the effective Higgs fields x and ^, on the branches parametrised by a. These 
display vortex ring structures described in detail in Section 3, similar to the vortex rings 
discovered previously [20j in the 3 + 1 dimensional Yang-Mills-Higgs system. 
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Figure 1. The mass M, the shooting parameter b and the value of the dilaton at the origin 
0(0) are shown as a function of a for n = 1, m = 1 spherically symmetric YMd solutions. (Here 
and in Figure 2, we use a normalization such that the mass of a = self-dual p = 1 YM solutions 
is M = q = n^.) 




Figure 2. The mass M and the value of the dilaton field at the origin (f)(0) are shown as a 
function of a for n = 1, m = 2 YMd solutions. 



18 



0.15 



0.1 



c^" 



-0.05 



-0.1 





e=o 




e=7t/8 




e= Tt/4 




e= 371/8 - 




e= 71/2 


m=l n=3 a= 0.21 - ' 

1 1 1 



-0.5 0.5 1 1.5 

logio(r) 

Figure 3a. 



c3 






1 


1 


1 








m=l n=3 oc= 0.21 


0.5 


9=0 














-1 


9= 71/8 








9=7t/4 






1.5 


9= 371/8 






-2 


9=7C/2 








1 





-0.5 0.5 1 1.5 

logio(i') 



Figure 3b. 
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Figure 3h. The YM gauge functions, the dilaton field and the topological charge density p 
are shown as a function of the radial coordinate r for a typical m = 1, n = 3 YMd solutions with 
a = 0.21. 
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m=2, n=2 upper vortex branch: |^| at a=0.15 




Figure 4a. 

m=2, n=2 lower 2-nodes branch: |^| at a=0.2 
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Figure 4b. 
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m=2, n=2 lower 2-nodes branch: Action density at a=0.2 
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m=2, n=2 upper vortex branch: Action density at a=0.2 

Figure 4c. 




Figure 4d. The modulus of the effective Higgs field ^ is shown for the uper branch m = 2,n = 2 
solutions at a = 0.15 (vortex, Figure 4a) and a = 0.20 lower energy branch solution (double 
node, Figure 4b) as functions of the coordinates p and a. The action density distributions of 
these m = 2,n = 2 solutions at a = 0.20 are also shown on the lower branch (Figure 4c) and on 
the upper branch (Figure 4d), respectively. 
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Figure 5a. 
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Figure 5b. The values of the dilaton function at the origin (Figure 5a) and the unrescaled 

mass Ma^ (Figure 5b) of the configuration with m = 2, n = 1, 2, 3 are shown as functions of the 
effective couphng constant a. The location of the vortices as given by (pojCTq) and the position 
of nodes (2^0^^ z^^) is also presented in Figure 5a for m = 2, n = 2 solutions. 
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Figure 6a. 



0.25 
0.2 
0.15 
0.1 
0.05 
c^" 
-0.05 
-0.1 
-0.15 
-0.2 
-0.25 



m=2n=3 a=0.17 



-0.5 



0.5 
logio(r) 



e= 





e= 


71/8 


e= 


71/4 


e= 


371/8 


e= 


nil 



1.5 




Figure 6b. 
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